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Abstract 



X> ■ 

The proposed yield criterion depends upon two material constants and is proven to be smooth and 
convex under a simple condition. These properties induce a mathematical robustness that allows a 
further use in a damage mechanics model. The analytical gradient and hessian are given. The obtained 
yield surface is relevant to Kupfer's biaxial testings on concrete. The identification procedure, with 
respect to the classical uniaxial tension and compression testings, is detailed. 



Jo Resume 

Un nouveau critere de limite d'elasticite pour les betons. Le critere propose utilise deux 
constantes materiaux. Son gradient et son hessien sont donnes, sa convexite est demontree, cette 
derniere ne dependant que de la positivite des constantes materiau. Ces proprietes induisent une 
•^h ' robustesse mathematique qui autorise son emploi futur au sein d'un modele d'endommagement. La 
J>y surface de charge obtenue est en bonne adequation avec les essais biaxiaux de Kupfer. La procedure 
^ d'identification, a partir des contraintes limites d'elasticite identifiers en traction et en compression, 
i ^~t " est donnee. 
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1 Introduction 



Kupfer P performed biaxial testings on concrete (<ti, an, am = in principal stresses) that are 
recognized as a reference and relevant for the stress states met in classical civil engineering. Although 
his best known results are the peak stress curves, the figure 15 of this article exhibits the elastic domain 
of concrete and the present work is referred to it. Kupfer measured the elastic limit from the loss of 
linearity of the stress to strain curve: this identification depends upon the precision of the measure. 
Some authors obtain, in pure tension, a yield stress very close to the peak stress i.e. a fragile behavior; 
on the contrary, Kupfer identified the yield stress in tension approximately as the half of the peak 
stress. A similar result is also obtained by Terrien [2] thanks to a very accurate method, confirming 
Kupfer's results. 

Most of old elasticity criteria, such as Mohr-Coulomb's one, or more recent approaches [3], are based 
on stress vector considerations. They may be accurate, have a strong mechanical sense and allow to 
define the crack orientation, but are non smooth. These corners are not experimental evidences and 
lead to numerical difficulties when the criterion is used in a damage model. 

Actually, the most used criteria are simple functions of the stress tensor a and the stress deviator <r d 
classical invariants I\ = tr(cr), J' 2 = tr(a d .a d )/2 and J3 = tr(cr d . )/3 whose polynomial structure 
insure the smoothness of the yield surface. Among them, we shall compare to the present model the 
De Vree jlj, Willam [5] and Ma'iolino [6] expressions (all using two constants). The identification will 
be made preferentially with respect to the uniaxial compression and tension yield stresses, respectively 
denoted by points C and T on figure (JTJ) - 

The De Vree criterion depends upon (Ji, J'^); its identification, with respect to the points C and T, 
gives the constants (7 = 7.75, k = 3.9 10 -5 ). The Mai'olino's criterion has been recently proposed for 
rock materials; it takes into account the effect of the third invariant J3 and, in order to keep a two- 
constants form, a linear dependancy {a + = I\ m — I\) has been retained here. The model, in this form, 
cannot fit both points C and T (the identified constant L s = 0, 5 reaches the bound allowed for the 
convexity of the criterion) then the identification has been done with respect to the points C and P, 
giving a good compromise and leading to the constants (ii m = 2.40 MPa, L s = 0,62). The Willam's 
criterion also makes use of the third invariant of the deviator J3. The identification with respect to 
the points C and T leads to the constant e = 0.5 which is the limit for the convexity of the model 
(dashed curve in figure [T]); again, a better identification has been found from points C and P (plain 
curve in figured]), with the constants (e = 0.68, m = 5.3 MPa). It can be seen that these models (and 
many others, see [7] for a more complete comparison) fail to describe precisely the particular shape of 
the concrete's yield surface. 

2 The proposed criterion 

In compression, i.e. when (every) a\ ^ 0, the yielding of concrete involves a diffuse microcracking, 
and high stresses. The mechanism in tension, i.e. when (at least one) a\ ^ 0, is very different, with low 
stress levels and a quick localization of the microcracks [8] . A part of the present criterion is related 
to the von Mises norm of the stress deviator that appears to be relevant to the compression states, in 
the low confinement range. The other part uses a stress tensor exponential (recalled in equation fl~2"I) 
in order to describe the dramatic effect of the positiveness of a principal stress, in the sense of the 
Rankine criterion, via an euclidean norm 

(I \A\ I = (AjAj) 1/2 )- The constant a y possesses a limiting 
role; the second constant o"o rules the dissymmetry of the yield surface and limits the numerical value 
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of the exponential term. The term V3 makes the second member equal to zero when a = 0. 
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Rewritten in term of principal stresses {a\, an, am), the equation ([2]) exhibits clearly the role of the 
second member that is positive if at least one principal stress is positive. One can remark that the 
numerical calculus of a tensor exponential does not require to diagonalize it (it is an integrated function 
in many softwares). 
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The constants (ao,a y ) can be identified with respect to the classical (and simplest) uniaxial com- 
pression and tension tests. Calling at > and a c < the experimental elastic limit stresses in pure 
tension and compression, corresponding respectively to points T and C on figure ([T]), the difference 
f(a t , 0, 0) - f(a c , 0, 0) = leads to: 



/(ff t ,0,0)-/(oc,0,0) = 
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CT0 V exp ) + 2 + "V I = ao v exp (^f ) + 2 ~ CTc v I (3) 

Each member is a monotonic function of <jq: their unique intersection, numerically calculated, gives 
the value of ctq. The equation f(a c , 0, 0) = 0, or /(at, 0, 0) = 0, gives the second constant a y . Here, the 
identified values are do = 0, 455 MPa and a y = 8, 00 MPa. The figure flU shows the good agreement of 
this criterion with the biaxial testings. The strong curvature around the point C of simple compression 
is well depicted. The figure P]), in the principal stresses space, shows that this criterion can be seen 
as a softened version of the combination of von Mises (for the cylindrical shape) and Rankine (for 
flattened faces) ones. The figure (J3j> shows the deviatoric sections of the criterion that evolve from the 
triangular shape of Rankine envelope to the circular shape of von Mises. 




Figure 2: Yield surface in principal stresses space 

(MPa). The plain lines correspond to the inter- Figure 3: Deviatoric view of the criterion, 

sections with planes 07 = 0. 



3 Gradient, hessian and convexity 

The use of a yield surface in a damage (or plasticity) model is helped by the knowledge of its 
gradient. We can already calculate the gradient of / with respect to the stress tensor a. The derivative 
of the first member of equation (pQ) is straightforward. For the second one, it is useful to use the 
property (|13p . The obtained gradient presents a simple (and intrinsic) expression: 

df_ g_ d , exp(2a/o- ) 
da \\a d \\ ||exp(a/ff )|| 

The hessian of / may also be useful in numerical calculus. For the first member of equation (pQ) we 
obtain the following expression, in which (g> represents the dyadic product and P d is the fourth rank 
projector onto the deviatoric subspace: 

a 2 lk d || f a d a d \ d 1 . j_ x x \ l xx ^ 

with P ijkl = — (dikOjl + Mjk) - o<Mkl (5) 



l^-Ml ~W\\ W\\J w ^ ll ^-2 KUikU ^ UilUik '~3 

The fourth rank operator between parenthesis in the equation © is a projector onto the four dimen- 
sions subspace of the second rank symmetric and deviatoric tensors orthogonal to a d ; its eigenvalues 
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are then (1,1,1,1,0,0) which proves the well known convexity of the von Mises criterion. For the 
second member of equation ([I]), it is convenient to have a change of variable a/ao = a': 
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Using the gradient (JU) leads to: 
d 2 \\exp(a!)\\ 
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Let us call iti the eigenvectors and A; the eigenvalues of a. Then iti are also the eigenvectors of a' 
and exp(a') (these tensors are coaxial with a) and \i/(Jq are the eigenvalues of a'. We consider now 
the orthonormal base of the second order symmetric tensors Bj_, with / G (1-.6), i.e. B_i : -Bj = <5n 
(the symbol ":" represents the tensor contraction (inner product) defined, for any symmetric second 
order tensors A and A', as A : A[ = AyA-: and 5 is the Kronecker delta). 
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it 2® it 2 , -85 = (lt 3 <g) itl + itl <g> ^3) 
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From formula (|14|) . the second derivative ([7]) expresses in the base Ri®Bj (in the Voigt notation, 
the terms 44, 55 et 66 would be divided by 2), in the general case (Ai 7^ A2 7^ A3), as: 
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From the convexity of the exponential function, the ry are positive if 00 is positive. The case A; = Aj 
is treated similarly from equation (I14p and it implies again ry > 0. The bilinear form associated to 
this expression can be written as follows, where A is a general symmetric second order tensor: 

S 3/2 A / " l|exp(a' 



. )(T r> , -A = Pi -A?i + P2A22 + PsA| 3 + 2r 23 A| 3 + 2r 3i A\ x + 2 n2 A 2 2 

+2p 2 p 3 (A| 2 - A 22 A 33 + Af 3 ) + 2p 3 p 1 (Al 3 - A 33 A n + A 2 X ) + 2p 1 p 2 (A 2 11 - A n A 22 + A^ 2 )(10) 

This expression is obviously positive and, added with the convexity of the von Mises term of the 
criterion, it proves the (strict) convexity of the proposed criterion /. This result can also be obtained 
without the hessian: it is easy to prove that the second partial derivatives d 2 f /da 2 are positive thus, 
from the Ball's theorem (number 5.1) [9], / is convex. 

4 Conclusions 



The proposed elasticity criterion offers a good description of Kupfer's biaxial testings on concrete. 
It depends of two constants that can be easily identified with respect to uniaxial compression and 
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tension curves. Its convexity is assured by a unique condition (o"o > 0), leading to mathematical 
robustness. Future work may consists in the development of a full damage model using this criterion 
as a yield surface. Another possible enhancement is to take the high confinement effects into account 
by replacing the von Mises part of the criterion by, for example, the Drucker-Prager or an elliptic 
closed expression. 



A Mathematics 



Let A = ^4ij"e^i <8> e j be a symmetric second order tensor expressed in an orthonormal basis {~<^i}; 
indexes (i,j,k, 1) refer to it. Let A p be the eigenvalues and {l^ p } the orthonormal basis of eigenvectors 
of A; indexes (p,q, r, s) refer to it. The derivative of the power low A™ expresses as: 
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else 

This allows to express the tensor exponential derivative as: 
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Two interesting properties are obtained from the (double) contraction of the previous expression of 
the exponential with A on the right side (i.e. A] 5 \d(exp(A))^/dA] i [) and with exp(A) on the left side 
{i.e. (eMA))ijd(eMA)h/dA kl ). 
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In the base of eigenvectors (T^ p <8> v q <8> u r ® ~^s)> the derivative (|12n has the following components: 
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if (p = q = r = s) 

if (p = r) / (q = s) et A p = A q 

if ( p = r ) / (q = s) et A p ^ A q 
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